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Abstract. Solar models computed with mass loss, micro- 
scopic diffusion of helium and heavy element, and with 
updated physics have been evolved from the pre-main se- 
quence to present day; they are compared to the observa- 
tional constraints including lithium depletion and to the 
seismic reference model of Basu et al. (1996), derived by 
inversion. Microscopic diffusion significantly improves the 
agreement with the observed solar frequencies and agree 
with the seismic reference model within ±0.2% for the 
sound velocity and ±1% for the density, but slightly wors- 
ens the neutrino problem. Neither microscopic diffusion 
nor overshooting explain the observed lithium depletion 
consistently with helioseismological constraints, while a 
mass loss process does it. Models computed with OPAL 
equation of state and opacities are in a better agreement 
with the seismic sound speed. To reach the level of pre- 
cision of helioseismological observations the accuracy of 
solar models still needs to be improved by one magnitude; 
any such improvement will necessitate equation of state 
and opacity data taking into account of detailed changes 
in the mixture. 
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1. Introduction 

The successful development of helioseismology imposes 
very strong constraints on the structure of the interior 
of solar models and has led solar modelers to improve 
the physical inputs of their models. One important step 
is the insertion of microscopic diffusion. The solar mod- 
els with microscopic diff usion (Proffitt 
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Brun et al. 1997 ) closely agree with constraints of the 
solar interior inferred by helioseismology as reviewed by 



Christensen-Dalsgaard et al. ( |1996| ). The recent improve- 
ments of equation of state (Rogers et al. 1996) and opac- 
ities (Iglesias & Rogers 1996[ ) have brought the difference 
between the solar sound speed of the seismic model of 
Basu et al. ( |1996[ ), from LowL data (hereafter, "seismic 
reference model") and the sound speeds computed from 
present day solar models within a rms discrepancy better 
than 0.2% (Cough et al. 1996); the values of the radius at 



the bottom of the convection zone (Christensen-Dalsgaard 
et al. 1991) and the helium abundance at the solar surface 
(Perez Hernandez & Christensen-Dalsgaard [1994 ; Antia & 
Basu 1994; Basu & Antia 1995) predicted by models agree 
within the error bar with their values inferred from helio- 



seismology (e.g., Basu 1997 ). Never theles s the large ex- 
cesses of neutrino rates (e.g., Bahcall 1997 ) and of lithium 
abundance predicted by the models still appear in strong 
conflict with the good agreement between the structure 
of the model and the inferred sound speed and density. 
Up to now, all attempts to mix the Sun in such a way 
as to minimize the discrepancies between solar observa- 
tions and predictions for, either neutrinos or lithium de- 
pletion, are ruled out by helioseismology. So mixing may 
be not the relevant process; it appears thus that for the 
neutrino problem the discrepancies are in pass to be un- 
derstood rather through a modest extension of the elec- 
troweak theory (Bahcall & Krastev 1996) than through 
improvements or change of physics. Promising processes 
for accounting for the observed lithium depletion are mass 
loss (B oothr oyd et al. |1992| ; Guzik & Cox |l995| ; Henyey & 
Ulrich 1995| ) and turbulent diffusion induced by rotation 
(Charbonnel et al. 1992, Chaboyer et al. 1995, Richard et 
19961 ) . 
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In this work, we have reconsidered the effect of micro- 
scopic diffusion and discussed the estimation of the heavy 
element content used in opacities calculations. Starting 
from chemically homogeneous pre-main sequence models 
we have computed calibrated solar models with and with- 
out microscopic diffusion. We have considered two differ- 
ent mass loss laws, undershooting and overshooting of con- 
vection zones, and also examined the effects of different 
equation of state and opacities. However turbulent diffu- 
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sion induced either by rotation or internal waves (Schatz- 
1992 , Montalban & Schatzman 1996| ) is ignored. 
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The models have been compared to the seismic refer- 
ence model through the quantities (c© — c mo d i)/ c o and 
(P© - /°modei)/p©, c modc i, c Q , p modc i and p e are, respec- 
tively, the sound speeds and the densities of the model and 
of the seismic reference model. The oscillation frequencies 
of the solar models have been compared to the GONG 
data observations. Likewise, the neutrino fluxes expected 
for the three experiments and the lithium depletion at 
present age have been computed and compared with the 
measurements. These solar models are updated and im- 
proved versions of the models discussed in Berthomieu et 
al. ( |1993| ). 

In Section)! global parameters and helioseismological 
constraints are briefly recalled for references and nota- 
tions. In Section|| the physics is described. Outlines of 
numerical techniques used so far are given in Section || 
Results and discussion are presented in Section^] and the 
conclusions in Section |[ Two appendices are devoted to 
detailed discussions of heavy element abundances which 
enter in equation of state and opacities. 



Table 1. Solar data at present day; the units are, 1Q M g for 
the mass Mq, 10 10 cm for the radius Rq, 10 33 ergs -1 for the 
luminosity Lq, Gyr for the age t©, dex ( 1 H=12) for the surface 
depletion of lithium Li©, SNU for $Ga and <3>ci, ev day -1 for 
3>Ka, AtHz for the frequency differences Suq2 and 5fi3. 



Mq 
Rq 
Lq 

tQ 

(Z/X)q 

Yq 

Rzc I Rq 
Li© 

$Ga 
$C1 
<&Ka 



1.9891 ±0.0004 
6.9599 ± 0.001 
3.846 ± 0.005 
4.52 ±0.4 
0.0245 x (1 ±0.1) 
0.246 ; 0.249 
0.713 ±0.003 

1.16 ±0.1 

69±S:S 

2.55 ±0.23 
0.29 ±0.02 
9.01 ±0.05 
15.90 ± 0.08 



Cohen & Taylor 
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fraction of heavy element to hydrogen and, iv) the solar 
mass. 



2. Observational constraints and calibration of so- 
lar models 

For reference, the solar global parameters, at present, and 
the constraints inferred from helioseismology are given 
m Table were the standard notations and the units 
are recalled in the caption. Recent observed neutrino 
fluxes at earth from the Gallium $Ga, Chlorine <I>ci and 
Kamiokande $x a experiments are reported. The frequency 
differences 8v 02 and Sviz which correspond to mean val- 
ues of recent observations from the ground based networks 
IRIS and BiSON and from the spatial experiments VIRGO 
and GOLF are given (see section 4.2). There Yq is the sur- 
face abundance per unit mass of helium, determined by 
Basu & Antia (1995) with respectively MHD and OPAL 
equation of state, and Rzc is the radius at the bottom of 
the convection zone. The solar age given in Table 1 corre- 
sponds to the start of main sequence evolution. We have 
adopted an evolutionary time of 4.5 Gyr for models start- 
ing from the zero-age main sequence and 4.55 Gyr for mod- 
els including pre-main sequence. The amount of heavy ele- 
ment per unit of mass for present day Z Q = 0.0181(1±0.1) 
is derived from Yq and from (Z/X)q the ratio of heavy 
element to hydrogen measured at the surface. 

The solar models discussed in this paper are calibrated 
within a relative accuracy better than 10 -4 by adjusting: 
i) the ratio I / H p of the mixing-length to the pressure scale 
height, ii) the protosolar mass fraction X p of hydrogen, iii) 
the protosolar mass fraction (Z/X) p of heavy element to 
hydrogen and, iv) the mass M p of the primitive Sun, in 
order that, at present day, the models have the values of 
Tablc[j]for, i) the luminosity, ii) the radius, iii) the mass 



3. Physics of Solar models 

3.1. Initial conditions and nuclear network 

Here the evolution of solar models, but one, include the 
pre-main sequence; they are initialized with chemically ho- 
mogeneous models, according to the procedure of Iben 



(1975) using a fixed contraction factor c = O.O^LqMq 1 K , 



see details in Morel ( |1997| ). 

The general nuclear network we used contains the fol- 
lowing twelve species 1 H, 2 H, 3 He, 4 He, 7 Li, 7 Be, 12 C, 
13 C, 14 N, 15 N, 16 and 17 which enter into the most 
important nuclear reactions of the PP±CNO cycles. The 
relevant nuclear reaction rates are taken from the tab- 
ulations of Caughlan and Fowler ( 1988] ). Weak screen- 
ing is assumed!!!. The abundance of any chemical is ex- 
plicitly computed i.e., no element is assumed to be a± 
equilibrium. The protosolar abundances of deuterium □ 
and of species heavier than helium are calculated us- 
ing their nuclide ratio s with respect to hydrogen (An- 
ders & Grevesse 1989 ); the protosolar mass fraction Y p 



of helium is taken as: Y p = 1 — X p (l ± (Z/X) p ); the 
protosolar 3 He abundance is derived from the protoso- 
lar 4 He = Yp abundance according to the isotopic ratio 
3 He/ 4 He = 1.42 x 1CP 4 . For lithium, we have used the 
ratio by number ( 7 Li/ 1 H) p = 3.277 dex; the unknown 
protosolar ratio by number 7 Be/ 4 H is taken to —3.58 dex 
i.e., about its equilibrium value at center of the Sun. As 



1 It has been recently demonstrated that weak screening is a 
very good approximate of the exact solution of the Schrodmger 



equation for the fundamental PP reaction (Bahcall et al. 1997). 
2 These computations do not use the re-estim ated value of 



the protosolar 2 H/ 1 H ratio by Gautier & Morel (1997) 
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a definition, we have considered that a pre-main sequence 
model becomes a main sequence model as soon as 99% of 
the energy released is from nuclear; this transition occurs 
at an age between 40Myr and 80Myr. At that time, due 
to the CNO nuclear energy produced by the conversion 
of protosolar 12 C onto 1 N, the models have a convective 
core of typical extent ~ 10% of the total radius and the 
outer solar convection zone has already receded about to 
its present day location. We emphasize the importance 
of pre-main sequence evolution for the study of the de- 
pletion of light element; with physical conditions of solar 
homogeneous zero-age main sequence and with the nuclide 
abundances of 2 H and 7 Li, the two reactions 2 H(p, 7) 3 He 
and 7 Li(p, a) 4 He are strongly far from equilibrium; a ho- 
mogeneous zero-age main sequence model can be found 
only if the abundances of 2 H and 7 Li are taken near to 
their equilibrium values, that is inappropriate for studies 
of light element depletion. 

3.2. The opacities and equation of state 

We have used the opacities of Livermore Library (Rogers 
& Iglesias 1992; Rogers et al. 1996) extended with the 
low temperature opacities of Kurucz's (1991) interpo- 



lated with a birational spline interpolation (Houdek & 
Rogl 1996) and, either the CEFF equation of state 
(Christensen-Dalsgaard & Dappen 1992| ), or the OPAL 
equation of state (Rogers et al. 1996] ). Because it depends 
not much on Z and to avoid too heavy numerical compu- 
tations, we have fixed the amount of heavy element in the 
equation of state to values closed to their mean amounts 
in calibrated models namely, Z = 0.0175 for solar mod- 
els without microscopic diffusion and Z = 0.0190 for so- 
lar models with microscopic diffusion. The opacities and 
equation of state are tabulated for mixtures with various 
amounts of hydrogen, Ai H , but with fixed ratios between 
the species of heavy element. During the evolution these 
ratios are modified by the thermonuclear reactions and 
by the microscopic diffusion and the new mass fraction of 
each isotope needs to be re-estimated as detailed in Ap- 
pendix^. In main sequence solar models a local maximum 
of 3 He occurs around R ~ O.3i?0, it amounts to 10% of the 
local helium content A3h c ~ 0.1 x A<i Hc ~ Z, here A3h c 
and A<i He are the mass fractions of, respectively, 3 He and 
4 He; as seen in Appendix^ this effect needs to be taken 
into account in the calculation of the amount of heavy ele- 
ment used for the determination of opacities and equation 
of state. 

3.3. Diffusion 

The turbulent diffusion is ignored in the solar models dis- 
cussed in this paper. Various treatment of the microscopic 
microscopic diffusion leads to dif fusion rates in good agree- 
ment (Bahcall & Pinsonneault 1995 ). We have used the 
microscopic diffusion coefficients of Michaud & ProfEtt 



( [1993 ). All chemicals including either A n cNO or Xz (see 
below), but 1 H and 4 He, are trace elements. Equation 
(14) of Michaud & Proffitt valid for a H-He mixture, is 
approximatively extended to a mixture including heavier 



elements, by writing (Alecian 1995) 



AiH^H = -(1 - - Z)ViYLc 



(1) 



instead of Jfi H Vi H = — (1 — Ai H )V4 He , therefore Ai H + 
A4 Ho + Z = 1 even with microscopic diffusion, here T^ijj 
and V4 Hc are respectively the diffusion velocities of hydro- 
gen and helium. 

The opacities and equation of state are functions of 
the heavy element content Z, through the number of free 
electrons and the abundances of efficient absorbers which 
do not necessarily belong to the nuclear network e.g., 56 Fe. 
Due to diffusion and nuclear reactions, Z changes as the 
Sun evolves, as well as the ratios between the abundances 
of chemicals; in the following we shall designate by Z K 
the amount of heavy element used for the computation of 
opacities in models and we discuss next three approximate 
estimates employed here. 

For the calculation of solar models without microscopic 
diffusion either, i) Z K is given by Eq. (A2) resulting of de- 
tailed calculations of abundances including the changes 
due to the nuclear reactions - as detailed in Appendix^ 
- but, the ratios between the abundances are not consis- 
tent with the values used for the creations of equation of 
state and opacities data or, ii) Z K is kept to its protoso- 
lar value Z K = Cte = Z p i.e., one neglects the changes of 
abundances of heavy element. 

For solar models with microscopic diffusion the alter- 
native is more subtle, in the core it is expected that the 
abundances of the heaviest chemicals are magnified by 
gravitational settling with the consequence of an enhance- 



ment of the opacity there (Turck-Chieze et al. 1993), there- 
fore a careful estimate of Z K is required to infer the effects 
of this process. 

A first possibility is to consider that, since opacity ta- 
bles are available only for one mixture of heavy element, 
we do not know how to take into account properly their 
changes due both to nuclear reactions and diffusion, so we 
just ignore them and take Z K = Cte = Z p . 

A second possibility is to model the diffusion of heavy 
element as a whole, by means of a mean fictitious chemi- 
cal Xz, with the mean atomic weight Mz, and the mean 
charge zz, of the heavy species; Xz being diffused as an 
extra trace element but it is not affected by the thermonu- 
clear reactions; then the ratios between the heavy element 
entering onto Xz are fixed, their changes due to thermonu- 
clear reactions are ignored and Xz can differ from its value 
given by Eq. (A2). 

We have also investigated a third possibility - labeled 
Z - namely, Z K is separated in two parts, the first consists 
of the chemicals heavier than helium which, belonging to 
the CNO nuclear network, are both diffused and nuclearly 
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processed, the second part consists of a fictitious mean 
non-CNO species X n cNO, of atomic weight M n cNO and 
charge z n CNO, which is only diffused (Eq. (A2)). Hence in 
the estimate of Z K , enter into account the changes of CNO 
abundances caused by diffusion, nuclear processes and the 
effects of the large gravitational settling of heaviest non- 
CNO species but, the ratios between the chemicals can 
differ from the values used for the computations of opac- 
ities and equation of state data. This third possibility is 
similar to the treatment of Proffitt ( 1994 ) , where in addi- 
tion, an attempt is made - namely its Eq. (3) - to approx- 
imate the changes of ratios between the heavy element in 
the calculation of Z K . 

For the mean chemicals, Xz and A n cNOi we have only 
taken into account the heavy element with abundances 
larger than 7.0 in the table 2 of Anders & Grevesse (1989), 
namely: C, N, O, Ne, Mg, Si, Fe; we have found, for X z : 
M z = 17, z z = 8 and, for A„ C no: -MnCNO = 27 and 

ZnCNO = 13. 

3.4- Convection, overshooting and undershooting 

The basic formulation of the standard mixing length the- 
ory is used. The convection zones are delimited according 
to the Schwarzschild's criterion. In convection zones and 
in their extensions by overshooting or undershooting, the 
chemicals are simultaneously homogenized and integrated 

199$ Morel 



with respect to time (Morel & Schatzman 



1997 ). The mixing includes every convective zone and its 



extent by undershooting or overshooting where the tem- 
perature gradient is kept equal to its adiabatic value (Zahn 



1991 ). The amount of overshooting of a convective core is 
taken as £ ov mm(H p , R cv ), Cov is the overshooting factor 
and i? cv is the radius of the convective core; at the outer 
limit i.e., in the atmosphere, no overshooting is allowed. 
With an overshooting factor of Cov ~ 0.3, the young Sun 
has a convective core, up to the age t ~ 1 Gyr; this results 
an enrichment in nuclear fuel which does not lead to sig- 
nificant changes in the present day solar model (Morel et 
al. |1996|). 



3.5. Lithium depletion and mass loss 

The depletion of lithium is mainly caused by the reaction 
7 Li(p, a) 4 He. For solar conditions (T ^ 1.5 x 10 7 K) it is, at 
least, 10 times more efficient than all the reactions involv- 
ing 7 Li, except 7 Li(d, n)2 4 He which is only two times less 
efficient (Caughlan & Fowler 1988), but the abundance 
of deuterium being almost zero, this reaction depletes no 
lithium. Therefore we can estimate that the abundance of 
7 Li resulting from our calculations is ~ 2% over-estimated 
- 6 Li destroyed around 2 x 10 6 K is not observed in the so- 



lar photosphere (Bohm-Vitense 1992) 



1993|) by the obser- 



As suggested (Soderblom et al. 
vations of the Pleiades (80-100 Myr) for stars with masses 
closer to the solar value only small lithium depletion seems 



to occur during the pre-main sequence while, from the 
observations of the Hyades (600 Myr) it appears that the 
most important amount of the depletion^ should occur 
during the first 500 Myr of the main sequence. Such 
a depletion can result from mass loss: in a star of mass 
M ~ I.IMq, at the onset of main sequence the convec- 
tion zone recedes at mass m > 1M©, the layers which 
correspond to the location of the present day convection 
zone i.e., m 6 [~ 0.97M©, ~ 1M©], are radiative and the 
temperature there is large enough to deplete the lithium. 
A mass loss of ~ 0.1M©, occurring at that time ejects 
the outer layers, and the present day convection zone is 
formed by mixing these lithium poor shells, with the result 
of the observed enhancement of lithium depletion. Follow- 
ing Guzik & Cox fll995p we have studied two laws of mass 
loss rate namely, a mild mass loss rate: 



M = -2 x 10~ 10 exp( 



t 



0.45 

and a strong mass loss rate: 



) M Q yr-\ 



M = -5 x 10~ 10 M© yr- 1 if M(t) > M©, 
M = otherwise; 



(2) 



(3) 



here t, in Gyr, is the age of the model and M(t) is its mass 
at time t; with the mild mass loss rate the protosolar mass 
M p ~ 1.1M© is adjusted to fit the present day solar mass. 



4. Computation of models and frequencies 

4-1. Solar models 

The solar models discussed in this paper have been com- 



puted using the code CESAM (Morel g_997j) ; al ong th e evo- 
lution atmospheres are restored (Morel et al. 1994 ). The 
number of mass shells and the time steps are fixed accord- 
ing to numerical criteria; at the onset of calculations (pre- 
main sequence) the number of integration points is typi- 
cally of the order of 550 and around 1030 for the present 
day solar models (~ 950 for the core and the envelope, 
~ 80 for the atmosphere); the whole evolution involved 
about 140 models, with 60 models from zero-age main 
sequence to present day. It has been checked that the nu- 
merical internal accuracy on sound speed is better than 
5 x 10~ 4 for the whole model, but differences below that 
level are perhaps meaningless due to uncertainties of in- 
terpolated physical data. It is assumed that the lost mass 
is simply detached from the star by some process which 
is not described. We also emphasize on the fact that in 
our calculations, the 7 Li depletion results from simulta- 
neous computations of chemical changes due to nuclear 
reactions, diffusion and convective mixing. The estimated 
fluxes for the three neutrino experiments, $ci, 'I'Ga and 
$Ka are computed according to Berthomieu et al. (1993). 



despite the higher metallicity of Hyades. 




Fig. 1. Fractional differences for sound velocities (left) and dens ities (right) between the solar models without microscopic 
diffusion of Table ^ and the seismic reference model of Basu et al. fll996| ). The main features of the physics used in each model 
are schematically recalled with the notations of Table|. SI {CEFF, OPAL92, PMS} (thin full), S2 {CEFF, OPAL95, PMS} 
(thick Ml), S3 {OPAL, OPAL95, PMS} (thin dashed), S4 {OPAL, OPAL95, ZAMS} (thick dashed). 



4-2. p-Mode Oscillation Calculations 



5. Results and discussion 



The frequencies of linear, adiabatic, global acoustic modes 
of the solar models have been computed for degrees £=0 
to 150 and have been compared to the observations. For 
the computation of frequencies the models are extended 
to about 1800 shells. The characteristic low degree p mode 
frequency differences Av n ,i = v n ,i — ^n-i.e+2 for £=0 and 
1, which provide information on the properties of the core 
of the Sun, have been fitted by linear regressions with 
respect to n: 



5v„ 



no), 



with no = 21 and for £ — 0,1, both for the observations 
and the theoretical frequencies. The quantities 8vo2 and 
(5^i3 given in Table 1 arc weighted mean values of the 8v n ^ 
der ived from BiSON (Chapli n et a l. |1997| ), IRIS (Gelly et 
al. [1997D, G OLF (Grec et al. |1997| ) and VIRGO (Frohlich 
et al. 1997). Despite the fact that solar activity is not ex- 



pected to have a large influence on the solar core structure, 
we have also represented in Fig. || another mean restricted 
to low activity observations. 

For the gravity modes which have not yet been ob- 
served, we give the characteristic asymptotic spacing pe- 
riod P according to Provost & Berthomieu (1986). 



We have computed solar models without microscopic dif- 
fusion, (labeled "Sn", n = 1, ...4) and solar models with 
microscopic diffusion (labeled "Dn", n = 1, . . . , 12) using 
the physics and physical data described Sect.|| and sum- 
marized in Table H for each model. 



5.1. The solar models without microscopic diffusion 

The global characteristics of solar models without micro- 
scopic diffusion are summarized in Table ||. With respect 
to our previous solar models without microscopic diffusion 



(Berthomieu et al. 1993) illustrated by SI, the changes are 
mainly due to equation of state and opacities. Figure]!] 
plots the comparisons with the seismic reference model; 
there the fractional differences for the sound speeds and 
densities are plotted with respect to the radius. Globally 
the models S are hardly closer than ~ ±1% for the sound 
speed and ~ ±10% for the density to the seismic refer- 
ence model. The comparisons between {SI, S2} and the 
seismic reference model, show that the improvement of 
models towards the seismic reference model is mainly due 
to the new opacities. 

Both frequency differences <5^o2 and 5v\3 are far from 
the observed values by more than 3tr - see Table|. The 
comparison with the GONG frequencies and the theoreti- 
cal frequencies for S3, plotted in Fig.[|, shows a noticeable 
dispersion of the curves for low and large degree modes, 
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Table 2. Summary of physical inputs for solar models with and without microscopic diffusion of Table g and Table |4| For each 
model the flag signals the physical option employed. In the first three lines, the evolutionary time is mentioned; the two next, 
indicate either, if the evolution includes the pre-main sequence and then the nuclear network with 12 species is used or, if the 
evolution is initialized from the homogeneous zero-age main sequence then it uses the simplified nuclear network with 2 H, r Li 
and 7 Be at equilibrium. In the four next lines the flags indicate if overshoot (Ov-shoot), undershoot (Un-shoot), mild mass loss 
rate (mMLR) or strong mass loss rate (sMLR) are allowed for. In the next three lines is indicated how the amount of heavy 
element is computed for the calculation of opacities (see text). In the last lines the mass ratio of heavy element to hydrogen, 
used for the calibration, the opacities and the equation of state used are given. 



SI S2 S3 S4 Dl D2 D3 D4 D5 D6 D7 D8 D9 D10 Dll D12 



age 4.65 Gy 
age 4.55 Gy 
age 4.50 Gy 
PMS 12 Chem. 
ZAMS 9 Chem. 
Ov-shoot 
Un-shoot 
mMLR 
sMLR 

Z K = Z P = Cte 

z K =z 

Z K = Xz 

(Z/X) e =0.0245 

(Z/X) 0=0.0260 

opa-OPAL92 

opa-OPAL95 

EOS-CEFF 

EOS-OPAL 



due to the bump (~ +0.6%) observed in sound speed be- 
low the convective zone on Fig.[l]. 



o 



_ o 



S3 


D10 


,"| . i v .'', 





3000 3000 4000 3000 3000 



Fig. 2. For models S3 and D10, the normalized frequency dif- 
ferences between the GONG frequencies and the theoretical 
frequencies are given for different degrees as a function of the 
frequency. For each degree, the points are linked by a contin- 
uous line: I = 2, 3, 4, 5, 10, 20 (dotted); £ = 30, 40, 50 (full) 
and £ = 70, 100, 120, 140 (dashed). 



Consequence of the pre-main sequence. In order to inves- 
tigate the effects of the pre-main sequence quasi-static 
contraction, the model S4 has been computed with the 
species 2 H, 7 Li and 7 Be at equilibrium, initialized as a 
homogeneous zero-age main sequence model and evolved 
4.5 Gyr i.e., ~ 50Myr less than the models including the 
pre-main sequence. As reported in Table || the global char- 
acteristics of S3 and S4 are similar and only a slight dif- 
ference on 5uq2 and 8v\^ reveals a very little change on 
the seismic properties of the models. The sound velocity 
of the two models are almost identical (see Fig.l), lead- 
ing to normalized differences of frequencies between the 
two models of the order of 0.5 /iHz. Similar comparisons^ 
performed for solar models with microscopic diffusion give 
differences of about 1 /iHz. Thus we can conclude that pre- 
main sequence evolution and general nuclear network are 
a necessity only for the study of specific relevant processes 
e.g., the lithium depletion. 



not detailed in Tabled for sake of briefness. 
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5.2. Standard solar models with microscopic diffusion 



Table 3. Global characteristics of solar models without micro- 
scopic diffusion. I is the mixing-length, Hp the pressure scale 
height; Y p , and X n cNO are the protosolar abundance in mass 
of respectively, helium and non-CNO heavy element; (Z/X) p 
is the protosolar ratio, in mass, of heavy element to hydrogen; 
LizAMS and Li0 respectively are the surface depletions in dex 
(H= 12) of 7 Li at zero-age main sequence and at present day; 
Yq, Zq and Rzc respectively are, at present day, the surface 
abundances, per unit of mass, of helium and of heavy element 
and the radius, in solar units at the bottom of the convection 
zone; T c , p c , Y c and Z c are the central values at present day re- 
spectively of, the temperature in units of 10 7 K, the density in 
g cm -3 the abundances, per unit of mass of helium and of heavy 
element. <3?Ga and $ci in SNU and 5>Ka in events day -1 , are 
the expected fluxes for the three neutrino experiments namely 
Gallium, Chlorine and Kamiokande; 81/02 and Si/13 are the av- 
erage values in piHz of the frequency differences between the 
radial p-modes of degree 1 = — 2 and £ = 1 — 3. Po is the 
characteristic spacing period of g modes in minutes. Data of 
Table|l] for the Sun (0) are also recalled. The physics used in 
the models is summarized in Table |[ 





© 


SI 


S2 


S3 


S4 


l/H p 




1.78 


1.80 


1.78 


1.78 


Yv = Y & 


0.246 


0.269 


0.271 


0.268 


0.268 


Z p = Zq 




0.0175 


0.0174 


0.0175 


0.0175 


(Z/X) p 




0.0245 


0.0245 


0.0245 


0.0245 


X n CNO 




0.0035 


0.0035 


0.0035 


0.0035 


LizAMS 




3.18 


3.15 


3.08 




Li© 


1.12 


3.18 


3.15 


3.08 




Rzc 


0.713 


0.733 


0.726 


0.723 


0.723 


T c 




1.535 


1.539 


1.538 


1.537 


Pc 




146.3 


147.1 


146.5 


146.3 


Yc 




0.615 


0.619 


0.617 


0.616 


z c 




0.0180 


0.0179 


0.0180 


0.0180 


$Ga 


79 


120 


121 


121 


121 


$C1 


2.55 


6.11 


6.30 


6.22 


6.20 


*Ka 


0.28 


0.47 


0.48 


0.48 


0.48 


8u 2 


8.99 


9.35 


9.30 


9.34 


9.30 




15.87 


16.42 


16.37 


16.40 


16.36 


Po 




36.51 


36.48 


36.45 


36.51 



The global characteristics of solar models with microscopic 
diffusion are summarized in Table || The comparisons of 
models with the seismic reference model are illustrated 
by Fig. H to Fig. ^ where the fractional differences for the 
sound speeds and densities are plotted with respect to the 
radius for the interval [0.05i? Q , O.85i? ] (Basu et al. 1996). 



The photospheric helium content of our models with 
OPAL equation of state is of the order of 0.245, smaller 
than the two determinations of Basu & Antia (1995). For 
models with the CEFF equation of state this helium con- 
tent is slightly increased. For all models D, except model 
D6 with penetration of convective elements in the inner 
radiative zone, the radius at the bottom of the convec- 
tion zone is within the error bars of the value inferred by 
helioseismology (see Table Q) in agreement with previous 
works. 

At the center the microscopic diffusion increases the 
density and reduces the available nuclear fuel, namely -"-H; 
a consequence of the calibration is a larger central tem- 
perature for solar models with microscopic diffusion than 
for solar models without microscopic diffusion; therefore, 
as noticed by numerous authors, the predicted neutrino 
rates which are mainly sensitive to nuclear parameters 
(e.g., Turck-Chieze & Lopez 1993; Dzitko et al. 1996) are 
slightly larger for solar models with microscopic diffusion 
than for solar models without microscopic diffusion by 
~ 7%, ~ 25% and ~ 25% respectively for $ G a, $01 and 

$Ka- 



n CD 



CD 

id L 



4^ 



SI O D5 □ 
S2 O DB « 
S3n A 



S4 * 
DlO 
D2 + 
D3* 
D4© 



DB 
D9 O 
DIOH 
Dll X 
D13-& 



9 9.2 
<5i/ 02 (/xHz) 



Fig. 3. Comparison of observed and theoretical estimations of 
the two mean frequency differences Svq2 and 5^13. 



In order to compare the seismic properties of the core 
of our models to those of the Sun, we have plotted on Fig.|| 
the values of the two mean frequency differences <5fg 2 and 
<5^i3. The mean observed values (see Table 1) and their 
corresponding errors are represented by the large box. The 
grey box corresponds to observations without IRIS values 
which have been obtained at larger solar activity. It is clear 
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that for our solar models without microscopic diffusion 
these quantities are too large compared to the observed 
values while models with microscopic diffusion are in bet- 
ter agreement with the observations. However, the models 
D5 (with mild mass loss rate) and D8 (with a larger age) 
have a too small 5^02- This quantity is larger than the 
observed one for models Dll and D12 where opacity is 
computed with a heavy element abundance Z^—Z. Both 
penetrative convection and overshooting of the core in the 
early stage of solar evolution do not modify these param- 
eters. Comparing the positions of models S2, D2 and Dl 
respectively to those of models S3, D3, D10, we note that 
the use of OPAL equation of state instead of CEFF equa- 
tion of state slightly increases the values of both 8v$2 and 
<5i^i3 . On the other hand, these latter models which have 
the best physics are in a better agreement with the seis- 
mic sound speed as will be seen later on. A larger age or a 
smaller value of Z/X could bring them in the observation 
box. 

Globally the characteristic period Pq for g modes are 
lower by about 1 minute for the solar models with micro- 
scopic diffusion compared to solar models without micro- 
scopic diffusion. 

As exhibited in Fig. 2 and in Fig.|] to ^, the mi- 
croscopic diffusion increases the agreement between so- 
lar models and the seismic reference model. The depen- 
dence of the normalized frequency differences between the 
GONG frequencies and the theoretical frequencies on the 
degree i is very weak, leading to a small dispersion of 
the curves in Fig. 2 (right). Globally the models D are 
close to the seismic reference model within ±0.2% for the 
sound speed and within ±1% for the dens ity. In agree- 
ment with previo us works (Gough et al. 199c , Christensen- 
Dalsgaard 1997 ) the sound speed is lower than in seis- 
mic reference model below the convection zone. On the 
contrary, all our solar models with microscopic diffusion, 
but D3 and D12, present a systematic minimum on sound 
speed differences around the radius R ~ O.5i?0. Moreover, 
below R ~ 0.4_Rq all our models except models Dll, D12 
have a sound velocity smaller than the seismic reference 
model. The negative minimum around R ~ 0.2i? Q which 
is present in most comparisons with inverse sound speeds, 
only appears as a small depression in our model D12 com- 
puted with OPAL equation of state and diffusion of Z with 
Z k = Z. 

Another consequence of the microscopic diffusion, as 
shown in Fig. [j], is that the evolutionary paths in the HR 
diagram of solar models with microscopic diffusion are 
shifted towards higher effective temperatures compared to 
the paths of solar models without microscopic diffusion. 

5.2.1. Sensitivity to equation of state and Z K 

By inspection of Figs 4, 5, 6, it is seen that the rela- 
tive sound speed differences with seismic reference model, 
have large variations of amplitude for models with CEFF 




3.B 



3.75 



3.7 3.65 

L °g T eff 



3.6 



3.55 



Fig. 7. HR diagram for models Dll {OPAL, sMLR, 2} 
{thin full), D12 {OPAL, 2} (thick Ml), Dl {CEFF, Z p } 
(thin dashed), D10 {OPAL, Z p } (thick dashed), S2 {CEFF, 
OPAL95, PMS} (thin dash-dot-dash), S3 {OPAL, OPAL95, 
PMS} (thick dash-dot-dash) and S4 {OPAL, OPAL95, ZAMS} 
(thin dotted); an enlargement of the main sequence is plotted 
in the panel. At zero-age main sequence the models have about 
the same luminosities, the effective temperatures for solar mod- 
els with microscopic diffusion exceed by ~ 2% those of solar 
models without microscopic diffusion. On the main sequence 
the evolutionary paths of S3 and S4 are superimposed. 



equation of state with a pronounced minimum around 
r = O.5i?0. The behavior is much smoother for mod- 
els with OPAL equation of state with a smaller mini- 
mum around r = 0.6Rq. This could be partly due to the 
way partial recombination in the solar radiative interior is 
avoided in CEFF equation of state. 

The model D10 (see Fig. obtained with the equation 
of state OPAL, is located closer to the seismic reference 
model than the model Dl obtained with the equation of 
state CEFF, but for the density profiles, Dl is significantly 
closer than D10 to the seismic reference model. The differ- 
ences in normalized frequencies are ~ 1/zHz between the 
two models. For the models D2 and D3 , the amount of 
heavy element Z K used for the computation of opacities 
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R / R^ R / R ( 



Fig. 5. Same as Fig.^ for the solar models with microscopic diffusion Dl {CEFF, Z p } (thin full), D6 {CEFF, Un-shoot, Z p } 
(thick Ml), D7 {CEFF, Ov-shoot, Z p } (thin dashed), D8 {4.65 Gy, CEFF, Z p } (thick dashed) , D9 {4.65 Gy, CEFF, Z p , 0.026} 
(thin dash-dot-dash). 




Fig. 6. Same as Fig.| for D3 {OPAL, X z } (thin full), D12 {OPAL, 2} (thick full), D10 {OPAL, Z p } (thin dashed), Dll 
{OPAL, sMLR, 2} (thick dashed) and Dl {CEFF, Z p } (thin dash-dot-dash); considering all constraints, Dll is our preferred 
model, but D12 without the lithium depletion constraint. 



in models is taken as the abundance of the mean fictitious 
chemical Xz- D2 is computed with the equation of state 
CEFF and D3 with OPAL. As seen Fig.g, for the sound 
speed D3 is closer than D2 of the seismic reference model, 
and also for density except for R & O.IRq. In agreement 
with Basu et al. (1996), however we find that the closest 
models to the seismic reference model are obtained with 
the opacities and equation of state OPAL. 

Comparisons of models D3 (Z K = Xz), D10 (Z K = 
Z p = Cte) and D12 (Z K = Z) allow a comparison between 
the three approximations used for Z K . As far as the sound 
speed is concerned, the models computed with an opacity 
value taking into account of the diffusion of heavy element 
are closer to the seismic reference model than models with 
Z K = Cte; the models D3 and D12 are very similar, but for 
R Si 0.2i? Q , D12 is closer to the seismic reference model. 
For the density, and R <; O.IRq the model D3 is closer to 
the seismic reference model than D10 and D12. Beneath 
the convection zone the sound velocity in model D12 is 
smaller than in D10 with a difference of the same amount, 
but of o pposi te sign, of those exhibited on Fig. 2(a) in 
Profntt (1994), however, the physical inputs of that work 
differ from ours. 

As a matter of conclusion, with the physical data used 
here, our study shows that the model computed with Z K = 
Z has a sound speed profile closer to the seismic reference 
model than the model computed with Z K = Xz, and that 
is reversed for the density profile. The model D10 com- 
puted with Z K = Cte, though farther from the seismic ref- 



erence model than D3 and D12, remains within ±0.2% for 
the sound speed and ±1% for the density. Therefore, using 
the seismic reference model, we are not able us to make 
a real discrimination between the three approaches em- 
ployed to compute Z K , nevertheless our "preferred" choice 
is Z K = Z. 

5.2.2. Sensitivity to age and (Z/X)q 

According to Table]!] the present day solar age is known 
within 0.4 Gyr, the model D8 has been evolved 0.1 Gyr 
more than the model Dl using the same physics; in the 
same way, model D9 has been calibrated for a larger value 
of (Z/X) e = 0.0260 instead of (Z/X) & = 0.0245 for 
model Dl. Figure^ shows that the sound speed in D9 is 
systematically smaller than in D8 by ~ 0.05% i.e., about 
our limit of accuracy, except around the center, while the 
density in D9 is greater in the convection zone and in the 
core. The spread in normalized frequencies is ~ 0.4/iHz 
between the two models. 

5.2.3. 7 Li depletion due to undershoot and mass loss 

All solar models with microscopic diffusion, except D4 and 
D5, reveal that the observed 7 Li depletion does not only 
result from microscopic diffusion, even if, as for model 
D6, an unders hootin g factor equal to the upper limit of 
Basu & Antia ( 1994 ), namely ( un = 0.1H p , is allowed for. 
The discontinuity of the temperature gradient due to un- 
dershooting is at the origin of spurious bumps, though 
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within the global accuracy, namely ±0.2% for the sound 
speed and ±0.5% for density, as exhibited in Fig.[|. But 
Ahrens et al. (1992) have shown from models without mi- 
croscopic diffusion that an undershooting factor £ un 0.3 
is needed to fit the present day observed lithium depletion; 
such a large amount is in conflict with the upper limit for 
Cim deriv ed from helioseismic observations ( Basu fc An- 
tia |1994 Monteiro et al. |1994| ; Provost et al. |l996|) . Our 
results for D6 (Table^) also show that the most impor- 
tant amount of lithium depletion occurs during the pre- 
main sequence. This is in conflict with the observations of 
lithium in Pleiades and Hyades (see Sect. 3.5). The lithium 
being mostly destroyed in the hot core at the end of the 
pre-main sequence (Morel et al. 1996| ) when the young Sun 
is still fully convective, a larger amount of undershooting 
will enhance the duration of this fully mixed stage and 
increase the pre-main sequence depletion. Therefore from 
the properties of our models it results that the 7 Li de- 
pletion observed at present day in the solar photosphere 
does not only result from overshooting and microscopic 
diffusion. 

In contrast in the models where the mass loss is in- 
cluded (models D4 and D5 of Table the lithium deple- 
tion occurs mainly during the main sequence. The compar- 
isons of theoretical sound speed and the density profiles 
with the seismic reference model (Fig.|J) show that the 
model D4 calculated with strong mass loss rate of Eq. (|J) , 
is closer to the seismic reference model than the model 
D5 computed with mild mass loss rate of Eq. (|J) . Indeed, 
stronger mass loss affects the structure of the present day 
model less because the mass loss phase is over faster. A 
model with strong mass loss rate becomes sensitive to the 
mass loss at the end of the pre-main sequence. There- 
fore, the evolutionary path in the HR diagram (Fig.]?]) of 
the model D4 reaches the zero-age main sequence with 
decreasing effective temperature and luminosity. Indeed 
both the strong and the mild mass loss rates look as rather 
ad-hoc assumptions so we have not attempted to adjust 
the free parameters of these laws in order to fit accurately 
the present day observed lithium depletion. Nevertheless 
it appears that a strong mass loss rate is a promising pro- 
cess to explain the observed lithium de pletio n in the Sun 
as already claimed by Guzik & Cox ( 1995 ) despite the 
large expected neutrino capture rates and the marginal 
agreement with the observed frequency difference <5fo2- 

5.2.4. Overshoot of the convective core of the young Sun 

The model D7 is computed with the physics of model Dl 
but with an overshooting of the convective core which ap- 
pears during the early stage of solar evolution by a factor 
Cov = 0.2. The global properties of these models (Tablc^) 
and the comparisons with the seismic reference model re- 
veal only slight insignificant differences at center between 
D7 and Dl (Fig.| and Fig.|) except a significant increase 
of the characteristic period Pq of the gravity modes. 



6. Conclusions 

We have computed solar models with our stellar evolu- 
tion code CESAM and using updated physics. Our mod- 
els are attempts to extend the study of the sensitivity 
of solar models with microscopic diffusion to pre-main se- 
quence, lithium depletion, mass loss, microscopic diffusion 
of heavy species, overshooting and undershooting. Effects 
of rotation and of turbulent diffusion are ignored here. We 
have compared the sound speed and the density profiles 
of the seismic reference model of Basu et al. ( 1996| ) to cal- 
ibrated solar model computed with various opacity data, 
equation of state, microscopic diffusion, mass loss rate, un- 
dershooting and overshooting amounts. The solar models 
with microscopic diffusion agree with the seismic model 
within ±0.2% for the sound speed and ±1% for density, 
while for the solar models without microscopic diffusion 
the agreement is hardly better than ±1% and ±10% re- 
spectively for the sound speed and density. For the solar 
models with microscopic diffusion the depth of the con- 
vection zone and the amount of helium at surface agree 
fairly well with their values inferred by helioseismology. A 
significant increase of the quality of solar models results 
from the recent improvements of opacities and equation of 
state, any amelioration of solar models will necessitate to 
take fully into account the changes of chemicals in opaci- 
ties and equation of state. 

Our models reveal that the lithium depletion observed 
at solar surface is certainly not due to undershooting at 
the bottom of the convection zone and can be explained 
by a strong mass loss occurring during the first 200 Myr; 
however, in the solar models described in this paper the 
turbulent diffusion induced either by the rotation or by 
the internal waves has not been taken into account. 

Considering all the observational and helioseismologic 
constraints our preferred model is Dll computed with 
strong mass loss rate, OPAL opacities and equation of 
state, microscopic diffusion of hydrogen, helium and, as 
trace elements, all the CNO species plus a fictitious mean 
no-CNO chemical which models the microscopic diffusion 
of the heaviest elements of the mixture; but if one ex- 
cludes the constraint on the lithium depletion, our pre- 
ferred model is D12 computed with the same physics but 
without mass loss. However, in the present state of art, 
the precision of solar models needs to be still improved by 
one magnitude in order to reach the accuracy of the solar 
data inferred by helioseismology. 
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Table 4. Global characteristics of solar models with microscopic diffusion. M p is the protosolar mass, £ 
are the overshooting and undershooting factors; the other captions are given Table|. 
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0.0191 
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0.0190 
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0.0191 


0.0194 


0.0196 


(Z/X) 


0.0271 


0.0271 


0.0271 


0.0270 


0.0267 
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0.0271 


0.0271 
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0.0270 


0.0276 
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0.0053 


0.0053 
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0.0052 
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0.0051 


0.0052 


0.0053 


0.0063 


0.0052 


0.0056 


0.0058 


LizAMS 
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3.13 
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3.00 
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0.247 


0.248 
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z e 


0.0179 
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0.0179 


0.0179 


0.0179 
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0.0180 


0.0189 
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2.92 
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2.17 
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Rzc 


0.710 


0.713 


0.711 
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1.560 


1.563 
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129 
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7.68 
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7.58 
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8.27 
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0.61 


0.62 


0.61 


0.66 


0.69 


0.61 


0.61 
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0.66 


0.60 


0.72 


0.66 
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9.00 


9.03 


9.10 


8.96 


8.85 


9.00 


9.00 
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9.04 


9.13 


9.09 


6V13 


15.96 


16.01 


16.04 


15.92 


15.77 


15.97 


15.99 


15.83 


15.92 


15.98 


15.98 


16.04 


Po 


35.73 


35.82 


35.87 


35.63 


35.27 


35.76 


35.97 


35.43 


35.65 


35.78 


35.76 
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A. Calculation of abundances 

The opacities and equation of state are tabulated for mix- 
tures with various amounts of hydrogen Xn , and heavy el- 
ement content Z; the ratios between the species of which Z 
is made of are fixed; along the solar evolution these ratios 
are modified by thermonuclear reactions and microscopic 
diffusion. Therefore, using the available data for opacities 
and equation of state, a difficulty is the estimate of Z as 
consistently as possible. At time t, Xi, the amount per 
unit of mass of the species labelled with i = 1, . . . , N c (N c 
is the total number of chemicals) is written: 



X t = 



and 



N e N 
1 = > , Xi = A" n CNO 



mM l m u 



X; 



(Al) 



here, rii is the number density, Mi is the atomic mass of 
Xi, Mi G 1R differs from the integer atomic number Ai, 
(Clayton 196S), m u is the atomic mass unit, N < N c is 
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the number of species entering into the nuclear network, 
-fnCNO is the amount, per mass unit, of a fictitious mean 
heavy element of mean atomic mass M n CNO which does 
not belong to the nuclear network and p is the density. 
Therefore Z is writen: 



Z = 1 — Xh — Xjjo — x 



for sake of clarity the labels "H" and "He" are used in 
place of integer indexes; Xr and Xn e also include, re- 
spectively, the isotopes of hydrogen and helium, (see Ap- 
pendix^). Following Fowler et al. (1975), we use the num- 
ber of moleg -1 , Yi = Xi/Mf, then the amount, per unit 
of mass, of the species indexed with i is written: 



X, 



YiMi 



YnCN0M n CN0 + EL Y 3 M i 



namely the abundances of hydrogen, and heavy element 
which are entries of opacities or equation of state, are re- 
spectively written: 



X H = 



Z = 



ynCNO-^nCNO 
ynCNO-^nCNO 



E, .• H it. ^ ; ^; 



Y, 



nCNO-WnCNO 



E 



N YM 

j=l 1 3 1V1 3 



(A2) 



therefore Eq. ( |A1| ) is fullfillcd, despite the fact that the 
quantity which is conserved is the nucleon number ex- 
pressed as: 



N 



i=l 



dt 



The similar approach of Richard et al. (1996) uses the 
normalization equation: 

1 = X u + Xt + Z. 



B. 3 He abundance, equation of state and opacities 

In main sequence solar models, a local maximum of 3 He, 
occurs around R ~ 0.3i?©, it amounts to 10% of the local 
helium content, there X3 Hc ~ Z; this effect needs to be 
taken into account in the calculation of Z K ; Fig. || exhibits 
the fractional differences on Z R , opacity and sound veloc- 
ity between the model SI used as a reference - for SI, 
Z K = Cte (see Tabled) - and solar models without mi- 
croscopic diffusion (not detailed here for sake of briefness) 
calculated with and without 3 He included into the helium 
content i.e., Z K = l- 1 H- 3 Hc- 4 He and Z K = l- 1 H- 4 He. 
From Fig.ra the fractional differences on sound speeds 
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Fig. 8. Relative differences for heavy element contents Z (top), 
k opacities (middle) and C sound velocities (bottom) along 
the normalized radius, between the model SI and solar models 
without microscopic diffusion in which 3 He either, is, Sy (full) 
or, is not, Sz (dash-dot-dash), included into the helium content 
for opacities calculations. 



are of the order of ±3 x 10 -3 i.e., of the order of accu- 
racy achieved with solar models with microscopic diffu- 
sion. Likewise Fig. ^| shows that the spread of O-C fre- 
quencies between sets of modes of constant I computed 
between GONG observations decreases by taking (left) 
Z = l- 1 H- 3 Hc- 4 He instead of (right) Z = l- 1 H- 4 He. 
The difference of the sound speed below the convection 
zone and around R ~ 0.3-R© seen in Fig.|| induces scaled 
frequency difference between the two models of the order 
of 1 /zHz, as seen in Fig. 10. 

Everywhere, but around ~ 0.3-R©, the small abun- 
dance of 3 He has no effect on opacities and equation of 
state; around R ~ 0.3-R© opacities and equation of state 
are mainly sensitive to free-free transitions, there 3 He and 
4 He are equivalent providers of electrons; therefore it is 
relevant to include 3 He in the helium content for the 
calculation of opacities. Similar phenomenon prevails for 
all components of Z which present noticeable changes of 
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Fig. 9. O-C frequencies between sets of modes of constant 
t computed between GONG observations and solar mod- 
els without microscopic diffusion computed respectively with, 
Z = l- 1 H- 4 He (Sz) and Z = l- 1 H- 3 He- 4 He (Sy). 
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Fig. 10. Scaled frequencies differences between the models Sy 
and Sz of Fig. ^| for modes of low degree (£= 0, 1, 2, 3, 4, 5, 10; 
dotted), intermediate degree (£= 20, 30, 40, 50; full) and high 
degree (£= 70 100 120 140; dashed). 

abundances due to nuclear reactions and diffusion e.g., for 
12 C and 14 N if R ^ 0.2i? Q . The sensitivity of the sound 
velocity to such effect being of same order of the accu- 
racy reached by solar models with microscopic diffusion, 
we emphasize that any improvement of that accuracy will 
necessitate opacities - even equation of state - taking into 
account the changes of mixture. 
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